A ring laser gyroscope without lock-in phenomenon 
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We theoretically and numerically study the effect of backscattering on rotating ring lasers by 
employing the Maxwell-Bloch equations. We show that frequency shifts due to the Sagnac effect in- 
corporating the effect of backscattering can be observed without lock-in phenomenon, if the strength 
of backscattering originating in the bumps of the refractive index is larger than a certain value. It 
is also shown that the experimental results corresponding to the theoretical ones can actually be 
obtained by using a semiconductor fiber-optic ring laser gyroscope. 

PACS numbers: 



I. INTRODUCTION 

Ring laser gyroscopes (RLGs) operate laser action 
on counter-propagating waves and measure their fre- 
quency difference (beat frequency) due to the Sagnac 
effect [ll y, Sj- The main problems for operation are 
caused by the effects of backscattering and nonlinearity 
inherent in a gain medium. For instance, when the ro- 
tation rate of the RLG is slow, frequency-locking effect 
induced by backscattering occurs between the counter- 
propagating waves and eliminates the beat frequency due 
to the Sagnac effect (Sagnac beat frequency), which is 
well-known as the lock- in phenomenon [3, y, [j, [f| Q- 
In conventional He-Ne RLGs, avoiding the lock-in phe- 
nomenon has been one of the important issues in order 
to measure the beat frequency at slow rotation rates, and 
many investigations have been intensively done for this 
purpose. 

While conventional RLGs utilize a gaseous medium 
(He-Ne) as a gain medium, a new type of RLG using 
solid states and a semiconductor material have recently 
received research attention 0, H ©, El EH El El El El 
llq |. As is well-known, the nonlinearity inherent in these 
gain media leads to a spatial hole burning effect and sup- 
presses one of the counter-propagating waves due to gain 
competition (in an ideal ring laser without backscattering 
and noise sources) , unlike in the case of gaseous medium 
0) Q • While several theoretical and experimental works 
exist for solidstate RLGs, many of the investigations have 
so far focused on overcoming or controlling such a non- 
linearity d EH El EH- Meanwhile, backscatter- induced 
phenomenon, such as lock-in, and the properties of the 
RLGs at slow rotation rates have not yet been examined 
in detail. 

In this paper, we theoretically and numerically study 
the effect of backscattering on the RLGs. In the descrip- 
tion, we use the Maxwell-Bloch equations, which consists 
of the Maxwell equations rewritten in a rotating frame of 
reference and the Bloch equations that have been used as 
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a simple model of solid-state gain media or as a semicon- 
ductor material. We show that Sagnac beat frequencies 
can be observed without the lock-in phenomenon, if the 
strength of backscattering originating in the bumps of 
the refractive index is larger than a certain value. In 
this case, we found that the rotation rate dependence 
of the beat frequency cannot be described by a conven- 
tional theory of the Sagnac effect, which gives a linear 
relation between the beat frequency and the cavity rota- 
tion rate [HQ- Instead, we demonstrate that it follows 
a novel formula of the Sagnac effect incorporating the 
effect of backscattering, which have been derived in Ref. 
[a, E3- I n addition, we report that such rotation rate 
dependence can actually be observed without the lock- 
in phenomenon in a semiconductor fiber-optic ring laser 
gyroscope (S-FOG) [HEo], which is composed of a semi- 
conductor optical amplifier as a gain medium and optical 
fibers to form a ring cavity. 

This paper is organized as follows: In Sec. HH we start 
with the Maxwell equations in rotating frames of refer- 
ence and review the eigenmodes of rotating ring cavities 
with backscattering to explain the Sagnac effect incorpo- 
rating the effect of backscattering. Then, the Maxwell 
equations coupled with the Bloch equations (Maxwell- 
Bloch equations) are introduced to describe the nonlinear 
interaction between the light field and a gain medium. 
We analyze the lock-in phenomenon by expanding the 
Maxwell-Bloch equations by the eigenmodes of a rotat- 
ing ring cavity with backscattering and derive the criti- 
cal backscattering strength, above which the lock-in phe- 
nomenon does not occur. In Sec. IIII1 numerical sim- 
ulations confirm whether the Sagnac beat frequency is 
observed without lock-in phenomenon. The experimen- 
tal results of S-FOG are reported in Sec. HVl Finally, the 
summary of this paper is provided in Sec. |V| 



II. THEORY 

Applying the electromagnetic equation of a naturally 
covariant form (or the general theory of relativity) to ro- 
tating ring cavities yields the Maxwell equations general- 
ized to a non-inertial frame of reference in uniform rota- 



tion with angular velocity vector [6|, [l8|, [lfj, |20| . Assum- 
ing that the vector is perpendicular to the plane where 
there is a ring cavity and that the light propagates one- 
dimensionally along the ring cavity, the following equa- 
tion can be obtained: 
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ds 
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2P(s)-E(s, t) + -— — P(s, t) + F!(s, t), (1) 



dt 



dt 2 ~ 



where s is the coordinate along the ring cavity in a ro- 
tating frame of reference with angular velocity $7. R(= 
L/(2tt)) denotes the radius of the ring cavity. P(s,t) is 
the polarization depending on the nonlinear response of a 
gain medium inside the cavity, whose expression is given 
in subsection III Bl Noise F\ is phenomenologically intro- 
duced in order to represent thermal fluctuation, quantum 
noise, cavity vibrations and so on; < F\(s,t) >= and 
< F 1 (s,t)F*(s',t') >= aiS(s ~ s')S(t - t'), where a x is 
the strength of the fluctuation, and < • • • > denotes a 
spatial-average or a time-average, c is the velocity of the 
light, and n(s) and /3(s) respectively denote the refrac- 
tive index and the background absorption. Throughout 
this paper, it is assumed that the effect of backscatter- 
ing from one of rotating waves to the opposite rotating 
waves can be described by the bumps of the refractive 
index and absorption. 



A. Eigenmodes of a rotating ring cavities with 
backscattering 

A key to understanding the effects of backscattering 
and rotation is to analyze the eigenmode structure of the 
ring cavities [171 |2l| . The eigenmodes in a ring cavity can 
be described by the following wave equation derived from 
the Maxwell Eq. JTJ) omitting absorption /3, polarization 
term P, and noise term F\: 



clockwise (CW-) rotating wave function and counter- 
clockwise (CCW-) rotating wave function with an iden- 
tical frequency ujq = (cfeo —)2nmc/(noL), where m is an 
integer (m = 1,2, ..). 

However, such modes with identical frequency cannot 
generally be described as eigenmodes when the spatial 
distribution of refractive index n(s) is non-uniform, i.e., 
there is a bump of the refractive index in the ring cav- 



ity: An (s)(= n 



7^ 0. In this case, the bump of 



the refractive index causes a linear coupling between the 
CW- and CCW- rotating wave modes due to backscatter- 
ing and the wavefunctions of the eigenmodes change into 
two standing waves [13, l2l l22| ■ Then, the degenerate fre- 
quency is split into two frequencies. For instance, when 
the bumps of refractive index An 2 <C 1 are so small, the 
frequency splitting Aloq[— Wi(0) — ^(0)] can be associ- 
ated with the strength of the backscattering which has a 
conservative nature, as follows fl7l.[21fl: 



A/o 



Aujq 
2vr 



tiottL 



V7 



(Hz), 



(3) 



where 7 is the intensity-reflection coefficient of the 
backscatters (7 <gC 1). Then the two standing wave modes 
with the splitting eigen frequencies are described as fol- 
lows; 



Ui(s) — VZcosnokos, U$ (s) = V2smnokos. 



(4) 



Moreover, when the cavity is rotated, the rotation af- 
fects the standing wave modes of a non-rotating cavity 
and causes shifts of their eigen frequencies. The effects 
can be shown by applying the perturbation theory for 
nearly-degenerate states typically used in quantum me- 
chanics to Eq. @ [23, [24], [23| . Based on the pertur- 
bation theory, the eigenmodes in a rotating cavity with 
backscattering can be represented as the superposition of 
the nearly-degenerate modes as 



U = c 1 U°(s) + c 2 U°{s) 1 



(5) 



if 
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„ RQw(Q) d TT . . n ,„, 
1 c 2 ds {S) = ,{) 



where E(s,i) oscillates with frequency u>(iY) as E(s,t) — 
U(s)e~ luJ ^ t + c.c and U is the wavefunction of the mode 
obtained under periodic boundary condition U(s + L) = 
U(s), where L is the cavity length. Below, it is as- 
sumed that the wavefunctions satisfy orthogonal relation 
1/(tiqL) <f UiU*n 2 ds = Sij up to the first order of Rfl/c. 
In the above, L is the total length of a ring cavity, and no 
is spatial-averaged refractive index n\ = 1/L § n 2 (s)ds. 

In the case of non-rotating ring cavity Q = 0, Eq. 
([2]) can be reduced to a usual wave equation (d 2 /ds 2 + 
n 2 (s)w J (0) 2 /c 2 )LT°(s) = 0, where uj(0) and U? respec- 
tively denote the eigen frequency and the wavefunction of 
a non-rotating cavity. In particular, when the refractive 
index is uniformly distributed inside cavity n(s) = no, 
it is known that the solutions of the wave equation are 



where the ratio of coefficients c\ and C2 is obtained as 



C2/C1 



sn 



(6) 



±\S 2 Q 2 + Aw 2 - Auj 



Then frequency luj(Q) of eigenmode j (j — 1,2) newly 
produced by cavity rotation up to the first order of R,Q/c 
is described as 



wi(fl) = uj + ^S 2 Q 2 + Au; 2 , 



iu 2 {n)=uj a -^S 2 Q 2 



Alu 2 , 



(7) 
(8) 



where Aloq is the frequency difference between modes of 
the non-rotating cavity and S is so-called scale factor, 
which can be represented as S = —?— \ § U® J^U^ds] « 

2ui R/(n c) = 4u A/(n cL) '2\. From Eqs. Q and 



((HI), we can finally obtain frequency difference A/n 
|Aw 2 (0) - Awi(n)|/(27r) as follows: 



A/n = ^ = vW + A/ 2 



(tfz), (9) 



where 5 = 5/(2tt) and A/ = Aw /(2tt). 

In conventional theory of the Sagnac effect, frequency 
difference proportional to the rotation rate has been de- 
rived between the CW- and CCW- rotating wave modes 
by assuming that there are no backscatters in a ring cav- 
ity pi 0- However, in the case of a ring cavity with 
backscattering, as seen in Eq. JjJJ), the relation between 
frequency difference A/n and rotation rate fl is not pro- 
portional, and it is modified by the effect of backscatter- 
ing 7. Then as seen in Eq. ((3J), the modes are expressed 
as superposed standing waves with the ratio shown in 
Eq. ([6]). This is the Sagnac effect in a ring cavity with 
backscattering. 



B. Model of a gain medium 

Next, we discuss the case of a ring laser taking into 
account the effect of a gain medium in order to examine 
how the Sagnac effect is observed in this case. 

To describe a gain medium in interaction with the light 
field, we use the well-known optical Bloch equations, 
which have been employed as a simple model of solid- 
state gain medium or as semiconductor gain medium: 



d_ 
dt 



p(s,t) = ~('j_L+iuj t )p(s,t) 

- iK,W(s,t)E(s,t) + F 2 (s,t), 



d_ 

dt 
2iKE(s,t)(p(s,t) 



W(s,t) = - 7// (W(s,t)-W 00 ) 

p*(s,t)) + F 3 (s,t), 



(10) 



(11) 



where it is assumed that the medium is not influenced 
by rotation. In the above, p and W denote microscopic 
polarization and population inversion, respectively. The 
two relaxation parameters, 7_j_ and 7//, are the transver- 
sal relaxation and longitudinal relaxation rates, respec- 
tively. u>t is the transition frequency between two level 
medium, and W^ is the pumping power. Noises F2 and 
-F3 are introduced phenomenologically in order to repre- 
sent thermal and pumping fluctuations: 



<Fi(s,t) > 

<F( Sl t)F;(s',t')> 



0, 

aiSijSis - s')S(t - t'), (12) 



where <Ji (i — 1,2,3) is fluctuation strength. The Bloch 
equations (fTT7|l and (jlip are coupled with the Maxwell Eq. 
(UJ) through the following relation between microscopic 
polarization p and macroscopic one P: 



P(s,t) = N(s)(p(s,t) + p*(s,t))Kh, 



(13) 



where N(s) is the atomic number density described as 
N(s) = N a Q(s), where 0(s) is a step function that is 1 
inside a region of a gain medium and zero outside it. 

Equations (JJ, (flu). (fTTj) . and $F5§ are the fundamen- 
tal ones describing the light field and the gain medium 
mentioned above. 



C. Two mode operation in rotating ring lasers 

In analyzing the dynamics of Eqs. fJJ, (fTO)) . (fTTj) . and 
(|13p in ring lasers, we expand these equations by the 
eigenmode basis of a rotating cavity derived in subsection 
III Al and derive the amplitude- and the phase- equations 
of the modes, although in conventional theoretical ap- 
proaches, CW- and CCW- rotating wave basis has been 
frequently used [2, Pi H , [ajp, HO, [HI [l2|, LL3| . As demon- 
strated in Ref. pj], [2l|, [22|, using the eigenmode basis 
makes it possible to simply discuss dynamical behaviors 
of ring lasers. 

To focus on the effect of a gain medium, we here omit 
noise terms F{ (i ~ 1,2,3) in the Maxwell-Bloch equa- 
tions fTJ), (fTO)) . (fTTj) . and (fT3")) , and assume that elec- 
tric field E(s,t) are described by two eigenmodes Uj 



(j = 1,2) as follows: 



c.c 



E(s,t) = E(s,t)e- luJot 



(14) 



+ c.c, 



where Ej and 4>j € Re denote the amplitude and the 
phase of the mode j of a rotating ring cavity, respectively. 
Then we assume that frequency ojq is close to the transi- 
tion frequency u>t of the gain medium, |o>o — ^t\/^t <C 1. 
In this case, solving Eqs. (JTOJ), (fTTj). and (fT3"|) 
by a perturbational method for small electric field 
yields P w -ia EW& + c.c and W = W«,(l - 
f3o\E\ 2 ) + 0{\E\ A ), where a = 2TrN a K 2 hu /j ± and 
00 = 4k 2 /(7^7//). By substituting them in the Maxwell 
equation (UJ and using the following orthogonal relation 
l/{nlL) § U*U jn 2 (s)ds = 5 lj + O (|i?fi/c| 2 ), where n is 
the spatial-averaged refractive index, n 2 , = \/L§n 2 ds, 
we can obtain the following equations for amplitude Ej ; 



d 

-r E i 
dt J 



(9j 



Sj E< 



C12E 2 



E, 



M c } cos * 



Mj sin ^ 



Mj-cos(2* + 6» x ),(15) 



where "net-gain" coefficient gj , self-saturation coefficient 
Sj, and cross-saturation coefficient en are given in Table 
UJ $ is the phase difference, i.e., $2 — <pi- Modulation 
amplitudes M c , M s , and Mj are described as: 

Ml = (I712I cos6» 7 + 3|^| cosOjE] 



-l^jlcosOz-jEi^E 



'3-JJ 



M; = (-l) 3 ^'(| 7l2 |sin( 



-le 



3-j 



> 1 + 3\£ J \sm9 J Ej 



)E,, 



(16) 



(17) 



Mi = \x\E 2 3 ^E 3 , (18) 



This condition is rewritten using Eq. ([3]) as follow: 



where 9y, 8j, and 9 X are the arguments of 712, £j, and \y 
respectively. 

On the other hand, the dynamics of phase difference 
*(= 4>2 — </>i) is described as: 



7>7c 



at 



(19) 



where Ao>n(= u>2(fl) — wi(f2)) is the frequency difference 
between the two modes, as shown in Eq. j9|. H, which 
is a periodic function of $ (period 27r), can be described 
as follows: 



H(9) = I712 



£ 2 + £ 2 2 
E1E2 

?2 , zt,2 



sin(\I> + 7 ) 



+ | X |(^ + £g)sin(2* + x ), (20) 

when amplitudes E\ and -E2 are small enough and £1 ~ 

Q « 1. 

The solutions of Eqs. (TT5)) and (fT5|) can be mainly clas- 
sified into two: one is a stable stationary solution with 
constant amplitudes and a constant phase difference, and 
the other is with time-dependent phase difference and 
amplitude modulating with the phase difference. The 
distinction can be clearly made by examining Eq. (|19p . 
which describes the time evolution of the phase differ- 
ence. 

Note that the form of Eq. (1191) resembles with the well- 
known locking equation [2J, [3|, 14, la, [6| ■ This means that 
the solution of Eq. (fT9|) markedly depends on the relative 
strengths of Aluq and max |if(\&)|, where max \H (W)| de- 
notes the maximal value of |.ff(\I/)|, which is obtained by 
solving the stationary equations of dEi/dt = 0. For in- 
stance, in the case of Aw^ < max|.ff(\IJ r )|, Eq. (|19p has 
stable fix points given by the solution of d^/dt = 0. 
According to definition of W, condition d^ jdt — im- 
plies that the frequency difference between modes 1 and 
2 equals zero, i.e., the lasing frequency of mode 1 becomes 
identical with that of mode 2, which is the lock-in (fre- 
quency locking) phenomenon. On the other hand, when 
condition Au>n 3> max|iJ(\I/)| is satisfied, the lock-in 
phenomenon does not occur. Thus, phase difference ^ 
can be approximately described as ^(t) « —Aunt + *0: 
where V^o is an initial phase, while the amplitudes Ej 
(j = 1,2) of the lasing modes oscillate with frequency 
Auin through the modulation terms of Eq. (fTS)) jlta ] . 

Here we emphasize that frequency difference Aujq in 
Eq. ([19")) not only depends on rotation rate £1 but also on 
the strength of conservative backscattering 7. As shown 
in Eq. §§§, Aluq never becomes zero even at rotation rate 
fl = when there are backscatters inside the ring cavity. 
This means that the lock-in phenomenon does not occur 
even at rotation rate fi = if frequency splitting Aloq is 
large enough to satisfy condition: 



•2.T-1 

Ac 2 



■Aw 



(22) 



Although criterion Alu c is too complicated to be analyt- 
ically expressed for general cases, it has the simple form 
for a case where the wavefunctions of the two modes can 
be expressed by Eq. (g]) because of 7 -C 1 and where 
wavelength 27r/fco of the light is shorter than the size of 
the ring cavity kg J Q(s)ds S> 1. Therefore g\ ~ gi = g 
and si ~ S2 = s. By solving equation dEi/dt = under 
such conditions, we can derive criterion Alo c : 



Aw r = 2 max 



\l\2\{K - l)sin* + #sin2* 



If + cos 2* 



(23) 



Aui > max|7J(*)| 



0=0 



Aoj r . 



(21) 



where K = s/x + 2 ~ 5. In the case of (dispersive cou- 
pling) 712 ^> g (linear net gain), since H($>) has max- 
imal value at \& m = ir/2(2m + 1), where m is an inte- 
ger, Eq. (|23|) is reduced to Au c = 2|7i 2 |, while in the 
case of 712 -C g, it is reduced to Aw c = 2g/^/K 2 — 1 
[26| . where H($) takes the maximal value at ty m — 
mix — 1/2 arctan \/K 2 — 1. 

So far, it has been believed that the lock-in phe- 
nomenon inevitably occurs at slow rotation rates in a ring 
laser with backscattering. However, this result means 
that, if conditions (f2"Tj) or (|2"2"|) are satisfied, the ring laser 
will have a property that the lock-in phenomenon can 
be naturally avoided without special techniques used in 
conventional RLGs, such as dithering 2]. It can be con- 
cluded that the frequency difference Awn described by 
Eq. © without the lock-in phenomenon can be observed 
even at arbitrary slow rotation rate. 



III. NUMERICAL SIMULATIONS 

In the previous section, we derived the amplitude and 
the phase-difference equations under two mode opera- 
tion by expanding the Maxwell-Bloch equations by the 
eigenmodes of a rotating ring cavity with backscatter- 
ing, and discuss the occurrence of lock-in phenomenon. 
Note that a perturbation theory for p and W and sev- 
eral approximations, including £, < 1, are employed in 
the derivation. Additionally, noise terms Fi are omitted. 
Thus, it is needed to check the validity of the theoret- 
ical predictions obtained in the previous subsection by 
the Maxwell-Bloch Eqs. flU, ([TO])- (fH)) including noise 
terms and full-order nonlinearity of a gain medium. In 
this section, we restrict ourselves to discuss the lock-in 
phenomenon by the numerical simulations. Detailed dis- 
cussions of the effect of noises on lasing dynamics will be 
reported elsewhere. 

Note that it will take a very long time to simulation 
the Maxwell-Bloch equations because the fast oscillation 
as well as the slowly varying envelop of the light field 
are simulated. The fast oscillation part of the light field 
has the constant frequency, which is very close to the 



TABLE I: Coefficients in Eqs. (15]), (16]), (17]), and 
(HU). 



Coefficient 



9 j = aotjWoo -in 

s ,^{e|[/,|V,i s 

2/3oVVoo rr^lrr |2irr |2„2j„ 
C12 = j? — J fc»|t/i I i- / 2 1 n as 

^ = ^-fe\U s \ 2 UtUS-jn 2 dB 

X = — 5"F- J 1 QUi Uin as 



«o = 



J "§£ 



2irN a K huo 

7-1 
-4t/ e|C/i| 2 n 2 ds 



"fij 



00 



■§[5U*Ujds 



4a n 



7x7// 
n 2 , = — § n 2 ds 



Physical Content 



Linear net gain 
Self-saturation 

Cross-saturation 
Modulation 
Modulation 
First-order factor 
Light-medium coupling 
Dispersive coupling 
Third-order factor 
Averaged refractive index 



transition frequency io t - Thus it is possible to separate 
the slowly varying part from the fast oscillation part in 
order to investigate the dynamics of rotating ring lasers. 
Let us suppose E and p to be the slowly varying en- 
velope of the electric field and the polarization field as 
E = Ec~ Wat + c.c and p — —ipe~ luJot , where the rota- 
tion wave approximation(RWA) is taken into consider- 
ation for the polarization. Applying the slowly varying 
approximation for time to Maxwell-Bloch Eqs. |T]), (flT)]) . 
(fTTj) . and (fTB"]) yields the following equations: 



8 F 



ds 2 



n 2 (s) 



2i- 



fln 8 



(-q no ds 

££ + »7(»)/5 + Fi(s,t), 



E 



di 



p = -(l± + i&o)p + kWE + F 2 {s, t), 



(24) 



(25) 



d_ 
Jit 



W = -1,,{W - Woo) - 2k{E~p* + c.c) + F 3 (s, i),(26) 



where space and time are made dimensionless by the 
scale transformation; (noUJo/c)s — » s and to^t — > t, re- 
spectively, and /3 = /3 c 2 /(n 2 .^), Cln = RVt/c, r)(s) = 
(2irN a Kh/nl)<d(s), jx = 7x/wo, 7// = 7//M), « = 
k/u>o, Aq = (u>o — Wt)/u>o- Fi is noise with the strength 
of fluctuation di = c 2 /(n^uj^ai. 

In our simulations, we chose a ring cavity with (dimen- 
sionless) radius R — 10 composed of a passive region and 
an active region where there is a gain medium. Then 
we suppose ri\ and ni to be the refractive index of the 



passive and active regions, respectively. Length I of the 
active region is set as I = (5V2 x 10~ 2 )L, where L is 
the total cavity length 27ri?. Then we set the values of 
most of the system's parameters as follows: k = 0.5, 
Woo = 2 x 10~ 3 , 7± = 5 x 10~ 2 , 7// = 10~ 3 , a = 10~ 7 
(j = 1,2,3), A = 0, and f3(s) = a L S(s), where 9(s) is 
a step function that is 1 inside the active region and zero 
outside it, and oil — 10~ 3 . 

Here, note that the conservative backscatters are in- 
troduced by the variations of the refractive index in this 
cavity. Since the strength of backscattering 7 is related 
to the variation as t/7 ex {n\/n\ — 1), the resonant fre- 
quencies of the eigenmodes of a non-rotating cavity are 
split into two, following Eq. (3]). The solid line in Fig. [1] 
shows the dependence of (dimensionless) frequency split- 
ting Awo/wo on the ratio of the refractive index n\/n 2 . 

For various values of n\/n\, we numerically integrate 
Eqs. ([24]) - (1 26 P to obtain two mode laser operation in the 
non-rotating ring laser (Q,o = 0). The detailed method 
of the numerical simulations is shown in Appendix [XJ 

In Fig. [1] we show the n\/n\ dependence of the fre- 
quency difference between the two lasing modes by closed 
circles. As seen in this figure, it turns out that there is a 
zone corresponding to the absence of the frequency dif- 
ference for n\/n\ < 1.006. This is because of frequency- 
locking between two modes. The regime (n 2 /n 2 < 1.006) 
corresponds to the case when strength 7 of conservative 
backscattering is smaller than critical value j c . On the 
other hand, the frequency difference can be numerically 
observed without the frequency-locking phenomenon for 
n\/n\ > 1.006, and its value approaches the solid line ob- 
tained from Eq (3]). Below, we refer to the two regimes, 
< n 2 /n 2 < 1.006 and n 2 /n\ > 1.006, as locking and 
unlocking regimes, respectively. 

Typical examples of the rotation rate dependence of 
the frequency difference are shown in Figs. Ufa) and (b) 
for the locking (n\/n\ — 1.005) and unlocking regimes 
{n\/n\ = 1.02), respectively. The frequency difference 
calculated by Eq. (9]) in each case is represented by the 
solid line. For the locking regime, when the ring laser 
is rotated but the rotation rate is small, the frequency- 
locking phenomenon is still caused and therefore, the fre- 
quency difference between the two lasing modes cannot 
be observed [Fig. [2] (a)]. When rotation rate SId is in- 
creased, however, the frequency difference start to in- 
crease and approaches the values calculated by Eq. (9]). 
The existence of the region corresponding to the absence 
of the frequency difference might be qualitatively similar 
with so-called deadband [2|. On the other hand, for the 
unlocking regime [Fig. [2] (b)], as mentioned above, the 
frequency-locking effect is not caused even at zero rota- 
tion rate. As the rotation rate is increased, the frequency 
difference following Eq. © increase without deadband, 
as shown in Fig. [2] (b). These numerical results confirm 
the theoretical prediction that frequency difference can 
be observed without the lock-in phenomenon in a ring 
laser with conservative backscattering. 
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FIG. 1: Frequency difference vs. ratio of refractive index 
712/ n\, where rii and ri2 are respectively the refractive index of 
passive and active regions in a non-rotating ring laser (Q — 0) . 
Frequency difference is made dimensionless as Auo/uio- Solid 
line denotes frequency difference between eigenmodes calcu- 
lated following Eq. © in the ring cavity (without a gain 
medium), while closed circles are results obtained by numer- 
ical simulations of Eqs. (|24p -(l26 p . Numerical results show 
that frequency-locking phenomenon occurs and frequency dif- 
ference becomes zero when n^/nl is less than ~ 1.006. 



IV. EXPERIMENT 

Lastly, let us examine if the theoretical results can ac- 
tually be obtained in a real experiment. For this purpose, 
we used a semiconductor fiber optic ring laser gyroscope 
(S-FOG) [la, LL6|. The experimental setup is illustrated 
in Fig. [3] The ring laser part of S-FOG consists of a 
semiconductor optical amplifier (SOA) as a gain medium 
and polarization-maintaining fibers to form a ring cavity. 
Conservative backscattering is mainly caused at the con- 
nection point between the SOA and the optical fibers in 
S-FOG, because there is an air region at the connection 
point, which leads to substantial changes of the refractive 
indices at the points [27j | . From the return loss measure- 
ments, the value of 7 was approximately estimated as 
10~ 5 < 7 < 10~ 4 . The extraction of the laser output 
signal from the ring laser part is achieved by a 95:5 cou- 
pler, where only 5% energy of the CW-output signal is 
extracted. It is detected by photo detector PD, as shown 
in Fig. [3J SPA denotes a rf spectrum analyzer that is 
used to observe the beat frequency signals of the output 
signal [28]. The entire system including power supplies 
and measurement instruments is mounted on a rotation 
table. 
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FIG. 2: (Dimensionless) rotation rate Qd vs. (Dimension- 
less) frequency difference Aajn/wo in rotating ring lasers with 
locking regime (nl/fti = 1.005) (a) and unlocking regime 
(nl/n? = 1.02) (b). Solid lines denote frequency difference 
calculated by Eq. <(9j in each case, while closed circles denote 
results obtained from numerical simulations of Eqs (|24|l 



A. Estimation of -y c 

To obtain the critical backscattering strength 7 C for the 
S-FOG, it is useful to estimate the parameters (dispersive 
coupling 712 and linear net gain g) that determine the j c - 
value (Eq. (|2^]l). 

Dispersive coupling 712 can be calculated on the 
basis of the definitions shown in Table Q] as 712 ~ 
T/(2nokl)\ sinnofcZsin6>| < T/(2nokl), where T is the to- 
tal loss (typically, 10 5 < T < 10 6 s~ 1 ) of the ring cavity, 
fco and I are respectively the wavenumber of the light 
and SOA length (I « 1.5 mm), and 9 is the parameter 
related to the position of the SOA. For the calculations, 
we used the wavefunctions (g]) because 7 <C 1. Using the 
S-FOG's parameters (lasing wavelength A ss 1563.1 nm, 
no = 1.445, and T ss 10 6 s _1 ) with the above equation 
yields dispersive coupling 712 < 10 4 s _1 . 

Linear net gain g of the SOA can be rewritten as 
g = go(a — 1), where a is the normalized pumping 
power defined by a = (W^ - W t h)fW t h (W t h is the 
pumping threshold), and coefficient go is estimated as 
90 ~ 1-3 x 10 7 (s~ 1 ) using the differential gain coefficient 
of the SOA. When we set normalized pumping strength 
a as a = 1.06 in the ring laser of fiber length L = 4.12m, 
linear gain coefficient g is estimated as g ~ 6.6x 10 5 (s _1 ). 
Using this value and values 712 = 10 4 s _1 , K = 5 
with Eq. (f2"2")l yields critical backscattering strength 
7 C « 7.2 x 10~ 6 , which is smaller than the measured value 
of 7(= 10~ 5 — 10~ 4 ). From these estimations, we expect 
that the S-FOG satisfies the condition (|2"2"|) and generates 
beat frequency signals without the lock-in phenomenon. 
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FIG. 3: Experimental setup of semiconductor fiber ring laser 
gyroscopes (S-FOG): SOA, a semiconductor optical amplifier 
(its temperature is kept at 25 degree by a thermo-controller) ; 
Coupler, 95:5 coupler ( 5 % of energy is extracted only from 
the CW-output signal inside ring cavity.); PD, photo detector; 
SPA, rf spectrum analyzer. 



B. Results and discussions 

Figure. [4] shows the measured rf spectra obtained 
from the PD output signal at normalized pumping power 
a = 1.06, where the spectrum is time-averaged to obtain 
clear beat signals. As predicted by the above estimation, 
a peak corresponding to the beat frequency can be seen in 
the rf-spectrum even at zero rotation rate (fi = Odeg/s). 
Then the peak frequency is changed as the rotation rate 
is increased. Fig. [5] shows the rotation rate dependence 
of the peak frequencies. The experimental result demon- 
strates highly linear characteristics only in cases of high 
rotation rate, while there is substantial deviation from an 
asymptotic line given by the dashed line for slow rotation 
rates. 

In addition, we confirmed that the similar rotation rate 
dependence with that shown in Fig. [5] can be observed 
even when changing the length of an optical fiber of the 
ring cavity part. According to Eq. ([3]), the beat fre- 
quency observed at zero rotation rate (fi = 0) can be 
controlled by changing the fiber length. In Fig. [6l we 
show the dependence of the beat frequency on the fiber 
length. We find that the dependence can be explained by 
using Eq. $5§ with the assumed backscattering strength 
7 = 7.0 x 10~ 5 , which is given by the solid line in Fig. [6] 
Assumed 7- value is consistent with the order of the value 
10~ 5 — 10~ 4 obtained from the return loss measurement. 
Moreover, to explain the rotation rate dependence of the 
beat frequency shown in Fig. [5] by using the assumed 7 
-value, we superimpose Eq. ^j with 7 = 7.0 x 10~ 5 in 
the solid curve of Fig. O It also turns out that the ten- 
dency of the experimental results is well reproduced by 
Eq. ([9]) . Consequently, we conclude that we observed the 
avoidance of the lock-in phenomenon for 7 > 7 C and the 
manifestation of the Sagnac effect © incorporating the 
effect of backscattering in the experiment of the S-FOG. 
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FIG. 4: rf spectrum at each rotation rate Q (deg/s). Sharp 
peaks, representing beat frequencies, can be observed in 1//- 
like noise. 
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FIG. 5: Beat frequency vs. rotation rate. Open circles denote 
beat frequency obtained from experiments, while dashed line 
denotes an asymptotic line. Eq. ((9| is represented by solid 
line. 



V. SUMMARY 

In summary, we studied the effect of backscattering on 
rotating ring lasers by using the Maxwell-Bloch Eqs. |T]), 
(fT0|) . (fTT|) . and (fl3|) and showed that the occurrence of the 
lock-in phenomenon can be characterized by the strength 
of the conservative backscattering that originates in the 
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FIG. 6: Beat frequency observed at rotation rate Q = 
(deg/s) vs. cavity length. Closed circles denote beat fre- 
quency at zero rotation rate. Eq. (|3} is represented by a solid 
line. 



bumps of the refractive index. Particularly, we revealed 
the existence of critical backscattering strength j c , above 
which the Sagnac frequency shifts described by Eq. © 
can be observed without the lock-in phenomenon. We 
also showed that the experimental results corresponding 
to the theoretical one is actually obtained in a S-FOG. 
We believe that the property that lock-in phenomenon 
does not occur is not unique for S-FOGs but univer- 
sal for RLGs using solidstate and semiconductor as gain 
medium. 
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APPENDIX A: NUMERICAL SIMULATION 
METHOD 



to carry out numerical simulations of dynamics. 

We divide Eq. (|24|) into wave propagation, polariza- 
tion, and noise parts. The time evolution operator of the 
wave propagation part is given as U T = e~ tHr , where H 
consists of the kinetic energy part T = —l/2(d 2 /ds 2 + 
2i$l]j/n Q d/ds) and the potential part including linear ab- 
sorption term V = — n 2 /(2n 2 t ) — if} as H = T + V. If time 
t is very small, the time-evolution operator can be split 
based on the symplectic-integrator method as follows: 



Ur 



_ ~iVr/2-iTT-iVr/2 



+ 0(r 3 ). 



(Al) 



The next time-step value of electronic field E(t + r, s) 
is calculated by operating the split U T on the previous 
time-step value E(t, s) and including the contribution of 
the polarization term and a noise term up to the first 
order correction of t: 

E(t + T,s)xi U T E(t, s) + T)(s)p(t, s)t + dfi(s, t), (A2) 

where dfx = Fit. 

Equation (|A2|) can be calculated by the follow- 
ing procedures; First, potential part e r lVT l 2 i s oper- 



ated on E(t,s), (Ei(t,s) 



-iVr 



E(t,s)). Then the 



Fourier translation of E\{t, s) is calculated as E\(t, k) = 
§ Exit, s)e- tks ds. Operating Ei(t, k) to e~ lf T yields the 
following result: 

E 2 (t,s) = — f e- l/2 ( k2+2nD/nok ) T E 1 {t,k)e lks dk.(A3) 
2vr J 

Then e _iW / 2 S 2 (t, s)(= E 3 (t,s)) is obtained. Finally, 
the next time step E(t + t, s) is obtained by adding po- 
larization term rj(s)p(t, s)t and noise term dfx(t,s) to 
2?3(£, s). dfi is generated by converting uniform random 
numbers generated by the M-sequence technique to nor- 
mal distribution random numbers by using Box-Muller's 
formula. 

For the Bloch equations, the next time-step value of 
polarization p(t + r, s) and population inversion Wit + 
r, s) are calculated by the Euler method: 

p(t + t, s) = p(t, s) + —p(t, s)t + df 2 {t, s), (A4) 
dt 



Let us explain the numerical simulation method of Eqs. 
Pi jl -pp ]) . The approximated Maxwell equation IJ21J has 
identical form to the Schrodinger equation with absorp- 
tion j3, polarization p, and noise iq; it allows us to use 
the techniques of the symplectic-integrator method [29| 



W(t + r, a) = W(t, s) + -W(t, s)r + df 3 (t, s), (A5) 
dt 

where dfi(t, s) — Fi(t, s)t (i = 2, 3) are also generated by 
the above method. 
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